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1 .  Introduction 

Many  magnetospheric  phenomena  of  the  Earth  are  directly  controlled  by  the 
solar  wind  which  flows  past  the  Earth.  Thus,  accurate  prediction  of  the 
interplanetary  conditions  of  the  Earth  is  a  prerequisite  for  a  workable  fore¬ 
cast  of  the  magnetospheric  dynamics  of  the  Earth. 

The  interplanetary  space  is  permeated  by  the  solar  wind  which  is  con¬ 
tinually  emitted  from  the  Sun.  Hence,  the  interplanetary  conditions  evolve 
as  the  coronal  conditions  undergo  changes.  Although  the  solar  wind  is  causally 
determined  by  what  is  present  at  the  solar  surface,  the  relationship  between 
the  interplanetary  conditions  of  the  Earth  and  the  coronal  conditions  of  the 
Sun  is  too  complicated  to  allow  a  workable  prediction  of  the  former  on  the 
basis  of  a  structureless  description  of  the  latter.  The  complication  results 
from  the  spatial  variations  as  well  as  the  temporal  evolutions  over  the  solar 
surface.  Evidently,  the  relationship  between  the  conditions  at  two  inter¬ 
planetary  points  linked  by  a  solar  wind  stream  is  much  simpler.  This  is  so, 
for  in  the  far  region  of  the  interplanetary  space  energy  deposition  by  means 
of  processes  other  than  thermal  conduction  is  insignificant.  Moreover,  most 
of  the  stream-stream  interactions  take  place  before  the  solar  wind  reaches 
the  far  region.  Therefore,  prediction  of  conditions  at  one  interplanetary 
point  on  the  basis  of  conditions  at  another  interplanetary  poir.t  v  more 
likely  to  be  workable. 

As  a  part  of  STIP  (Studies  of  Travelling  Interplanetary  Phenomena) 
program,  we  undertake  a  two-point  correlation  study  to  explore  the  feasibility 
of  a  prediction  scheme.  First,  we  elucidate  the  theoretical  solutions  that 
represent  the  solar  wind  streams  in  the  far  region  of  the  interplanetary 

space.  Then,  we  present  a  useful  approximation  for  the  explicit  expressions 

! 

of  these  solutions.  Such  an  approximation  provides  an  efficient  algorithm 
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for  routine  calculation  of  the  solar  wind  streams.  Finally,  we  use  the 
observational  data  from  the  Venus-bound  space  probe  Pioneer  12  and  the  Earth- 
bound  satellite  IMP  8,  in  conjunction  with  the  solutions  from  theoretical 
modelling,  to  discern  the  two-point  correlation  in  the  interplanetary  condi¬ 
tions.  Such  an  assessment,  aimed  to  ascertain  the  feasibility  of  the 
prediction  of  the  Earth's  interplanetary  conditions,  is  possible  with  the 
recent  availability  of  simultaneous  data  at  0.7  All  (where  Pioneer  12  was 
inserted  into  Venus'  orbit  in  December  1978)  and  1.0  AU  (where  IMP  8  was 
launched  into  a  geocentric  orbit  in  October  1973). 
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2.  Governing  Equations  of  Solar  Wind  Streams 

As  a  starting  point,  we  shall  write  the  magnetohydrodynamic  equations  in 

a  rotating  frame  of  reference.  Referring  to  a  heliocentered  spherical  coordinates 

(r,  8,  <j>)  which  rotates  at  a  constant  angular  velocity  ^  =  t  Q  ,  the  temporal 

z  o 

and  spatial  variations  of  the  physical  quantities  (viz,  mass  density  p, 
temperature  T,  thermal  pressure  p,  heating  flux  ($,  flow  velocity  u,  magnetic 
field  §,  electric  field  j£,  and  current  density  3)  in  rationalized  units  are 
governed  by  the  following  equations: 


■  +  •*  « 

,,U  +  7‘OU  =  0  , 

c  u 


(1) 


+  (u-v)ii  +  2i7xu  -  =  -Vp  +  pVy®+  J*  B  ,  ^ 


>t  °'  2^+  Ira?  - 


+  V .  o  ^u2+  -  — 3 +  pu  +  p^xS+Qi  =  0  ,  (3) 


+  VxE  =  0  , 


(A) 


V  *  B  =  0 


(5) 
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in  conjunction  with  the  following  relational  equations 


p  = 

K  pi  , 

(6) 

J  = 

(7) 

E  = 

-  UK  I  , 

(8) 

3- 

(  -K  vt)  •  g  _ 

- r -  B  • 

(9) 

g^ 

The  vector  q  is  the  axial  distance  from  the  solar  axis  of  rotation  (viz. 
q  =  r  sin  0).  The  constants  G,  M®,  K  and  p  are  the  gravitational  constant, 
the  solar  mass,  the  gas  constant  (viz,  Boltzmann  constant  divided  by  average 
mass  of  particles)  and  the  magnetic  permeability.  The  thermal  conductivity 
K  =  C  T5^2  depends  on  the  temperature,  with  the  coefficient  C  essentially  constant. 

We  remark  that  the  thermodynamic  quantities  p,  T,  p,  Q  as  well  as 
the  magnetic  field  B  and  current  density  J  are  invariant  in  the  coordinate 
transformation.  In  the  rest  frame  the  flow  velocity  is  u  +  and  the 

electric  field  is  E  -  fi^qlj,  x  hence  the  kinetic  energy  is 

~u2  +  +  ^Q2q2  and  the  Poynting  flux  is 

p  *E  x  B  +  (n^qp^B2)  1^  -  (fi^qp -1B(ji)B.  The  latter  is  clear  if  we  recall  that 
in  a  perfectly  conducting  medium  the  Poynting  flux 


E  *  B 


u  +  l 


■xr  'ii 


-  BB 


(10) 


consists  of  the  convected  flux  of  magnetic  energy  and  the  work  done  by  the 
magnetic  stress  (11  stands  for  the  unit  dyadic).  Moreover,  although  spatial 
derivatives  are  invariant  in  the  coordinate  transformation,  the  temporal 
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* 


t 


derivative  in  the  rest  frame  is  3/3t  -  QJD/Bi J>,  hence  an  azimuthal  variation  in 
the  rotating  frame  will  appear  as  a  temporal  variation  in  the  rest  frame. 


It  is  instructive  to  write  the  equation  of  azimuthal  motion  in  the  conser¬ 
vation  form  for  angular  momentum 


a  t 

ici  ?  s 


(11) 


P  q  P’-ijq  o  u,,  +  o  q':_P0  U  I 

=  -,4  +  “U_1B':)+  r-’ 

VP.'  .  / 

The  term  3/3t(pq20  )  +  V*(upq2Q  )  accounts  for  the  torque  2pQ  q  *  u  exerted  by 
the  Coriolis  force  and  the  term  3/3ij>  p  accounts  for  the  torque  exerted  by  the 
thermal  pressure.  Clearly,  qp(U.  +  qft  )  is  the  material  angular  momentum  in 
the  rest  frame,  and  (q-p -1B2) -  (p_1qB^)B  represents  the  flux  of  magnetic 
angular  momentum.  Also,  it  is  instructive  to  write  the  equation  of  energy 
conservation  in  the  form 


02) 


which  expresses  the  increase  of  entropy  S  =  K  log(T3/2/p)  by  the  heating 
fltiy. . 


b 


Now  we  consider  a  solar  wind  stream  that  appears  time-independent  and 
hydromagnetically  aligned  (i.e.,  the  flow  velocity  is  parallel  to  the  magnetic 
field)  in  the  rotating  frame.  Since,  by  equation  (8),  the  electric  field  is 
zero,  equation  (4)  is  satisfied.  Let  £  be  a  length  variable  that  increases 
outward  (viz.  away  from  the  sun)  along  the  axis  of  a  flux  tube.  Hence,  2  =  ±TgB 
(with  the  positive  sign  if  the  field  line's  polarity  is  outward  and  the  negative 

-V 

sign  if  inward),  u  =  l^u  and  Q  =  Equation  (9)  means  that  the  heat  flux 

conducting  along  the  flux  tube  has  the  magnitude 


-Cl 


5/2  dT 

6) 


(13) 


By  virtue  of  equation  (5),  the  field  strength  B  varies  inversely  with  the 
cross-sectional  area  of  the  flux  tube.  Accordingly,  equation  (1)  yields  the 
invariance  of  the  mass  efflux 


m  = 


and  equation  (3)  yields  the  invariance  of  the  energy  efflux 

I1‘  5  GM0  i  '  2  \ 

a  =  pj  V+  "<r  )  u  +  q 

C  "  B . 


(14) 


(15) 


along  the  flux  tube.  Furthermore,  if  the  azimuthal  gradient  of  the  hydromag- 
netic  pressure  is  negligible,  then  equation  (11)  yields  the  invariance  of  the 
hydromagnetic  angular  momentum  efflux 
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^  .  p  q  K +  gq) u  -  f1*  V 

B 


along  the  flux  tube. 


(16) 


In  our  exposition  of  solar  wind  streams  along  flux  tubes  that  spiral  on 
helioapexed  poloidal  cones,  the  equation  of  poloidal  motion  will  not  be  used. 
The  equation  of  radial  motion  can  be  written  as 


r  dur 
urdr 


\  P  r  dr*' 

_(yy  sln  9)  +  **<2  -  t  a? 


2  -TT' 


-  +  JL  J  R, 

r _ p  91 


(17) 


in  which  Jq  =  -p"1r-1d(rB(j))/dr  is  the  poloidal  current  if  the  azimuthal  gradient 
of  is  negligible.  The  equation  of  longitudinal  motion 


/ 

?iu 


du 

ar 


G.l© 

r 


(18) 


can  be  combined  with  equation  (15)  to  yield  the  entropy  equation 


4..  3kt  +  p  -1 A 

d  2  y  dj  o 


_  d  ft. 
d,'  ou 


(19) 


which  can  also  be  obtained  directly  from  equation  (12). 
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3-  Solutions  for  Spiraling  Solar  Wind  Streams 

In  the  far  region  of  the  interplanetary  space,  the  flux  tubes  are  likely 
to  spiral  on  poloidal  cones  (viz,  0  =  const.)-  This  means  Bg  =  0.  On  the 
other  hand,  the  magnetic  field  has  an  azimuthal  component 


B,  = 

? 


_  X  -TflCj  rasin'" 


r  sin  9 


(20) 


? 


-  l 


This  follows  from  equation  (16)  in  conjunction  with  the  kinematic  relationship 
Ug  =  Bg  u/B  implied  by  the  alignment  of  the  stream  tube  with  the  flux  tube. 

In  the  absence  of  stream-stream  interaction,  the  cross-sectional  area  of  a 
flux  tube  will  have  its  radial  projection  varying  as  r2.  Accordingly,  the 
invariance  of  the  magnetic  flux  along  a  flux  tube  takes  the  form 


y  -  r  B  cos  "X 


(21) 


in  which  the  spiral  angle  of  the  flux  tube 


B 

X  =  arc  tan  —  (22) 

B 


is  the  angle  between  its  axis  and  the  heliocentric  radial  line.  Next,  in  view 
of  the  kinematic  relationship  u^  =  Br  u/B,  it  is  seen  that  r2pur  is  equal  to 
'/,  '2.  Upon  the  use  of  dr/d2  =  cos  x  and  P  =<^?5yr2nr,  equation  (19)  can  be 


written  as 


9 


d 

"dr 


(  C  cosSt„d  t5/^  dT ^  3  dT  tn  ,  „  N  T  _ 

uf;  -  a if  ■ <2  +  x)t- 


(23) 


in  which  \  =  (r/u^) (du^/dr)  is  given  by  equation  (17). 


Now  specification  of  p,  u,  T,  B,  and  X  at  an  interplanetary  point  with 
radial  distance  r  will  determine  the  geometry  of  a  spiraling  flux  tube  passing 
through  that  point  as  well  as  the  physical  profile  of  the  solar  wind  stream  in 
that  flux  tube.  The  heat  flux  must  take  on  a  certain  compatible  value  deter¬ 
mined  by  the  above-mentioned  specified  values,  so  that  the  pressure  will 
vanish  at  infinity.  Equations  (13)  and  (23)  constitute  a  second-order  system 
of  differential  equations  for  T  and  Q.  It  determines  two  functional  combina¬ 
tions  which  are  invariant  along  the  flux  tube.  One  of  these  two  invariants 
must  take  on  a  certain  value  (say,  zero),  thus  providing  the  compatibility 
constraint  which  will  determine  the  heat  flux.  The  other  invariant  can  take 
on  an  arbitrary  value,  say,  ^ .  We  do  not  know  the  mathematical  expressions 
nor  the  physical  meanings  of  these  two  obscure  invariants.  At  any  rate,  the 
six  invariants  (namely ,“57.,  C,«C,  7,  J,  and  zero)  will  determine  the  six  varia- 
b1  '  (namely,  p,  u,  T,  B,  X  and  Q)  as  functions  of  r.  The  difficulty  lies  in 
the  judicious  determination  of  the  compatible  heat  flux  in  terms  of  the  observed 
values  of  the  remaining  five  quantities.  Since  this  is  a  two-point  boundary 
value  problem,  trial-and-error  method  of  numerical  shooting  in  conjunction 
with  numerical  integration  of  differential  equations  is  rather  tedious.  In 
the  next  section  we  shall  obtain  analytic  expressions  which  are  satisfactory 
approximations  for  the  compatibility  constraint  and  the  invariant 
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U.  An  Approximate  Solution  in  Explicit  Expressions 

We  have  formulated  a  deterministic  set  of  equations,  whose  solutions, 
satisfying  the  condition  of  vanishing  pressure  at  infinity,  will  describe 
solar  wind  streams  in  the  interplanetary  space.  The  solutions  of  the  differen¬ 
tial  equations  (13)  and  (23)  can  be  visualized  as  curves  in  the  (r,  T,  Q) 
space.  Not  all  solutions  are  physically  admissible.  The  admissible  solution 
should  have  the  property  that  T  ■*  0  as  r  •*  Some  properties  of  the  admissible 
solutions  are  clear  from  the  invariants  along  a  flux  tube.  First,  the  radial 
velocity  attains  a  finite  value  at  large  heliocentric  distances.  This  follows 
from  the  fact  that  its  kinetic  energy  is  converted  from  the  finite  amount  of 
thermal  energy  and  heating  flux.  Thus,  the  dimensionless  ratio  \  diminishes 
to  zero.  In  fact,  most  of  the  radial  acceleration  takes  place  during  the 
transonic  transition  which  occurs  in  the  corona  before  the  solar  wind  reaches 
the  interplanetary  space.  Secondly,  the  spiral  angle  of  the  flux  tube 
approaches  90°  at  infinity.  This  follows  from  the  variations  that  the  radial 
field  decreases  like  r  ~2  whereas  the  azimuthal  field  decreases  only  like  r-1. 
Indeed,  x  is  rather  large  and  increases  very  slowly  in  the  interplanetary 
space.  This  can  be  seen  from  the  equation 

tan  r  W  sin  "6  -  <£  cos  -  J 

tan  \  '  7oWn.rf^  -TfWvi  cos  X  -J  ’  <24> 

which  indicates  x  ~  actanO  r/rQ)tan  XQ)  in  the  interplanetary  space  where  the 
radial  velocity  does  not  change  much.  We  use  subscript  "o"  for  the  interplan¬ 
etary  "base"  at  which  data  are  specified. 
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It  is  reasonable  to  analyze  the  problem  by  approximating  the  slowly 
varying  X  and  X  in  equation  (23)  with  constants.  The  approximating  values  X 
and  x  may  be  chosen  as  and  XQ  at  the  specified  interplanetary  point.  In 
terms  of  the  three  similarity  variables  (Yeh,  1971) 


R  =  C  r  , 


C  'X 

X  = - — r  T' 

K  W7 


Y  = 


]_  C"'H  V 

k  in  j-  _"“ 


(25) 

(26) 

(27) 


the  approximated  equation  (23)  can  be  split  into  two  coupled  first-order 
differential  equations 


dR  _  dX  _  dY 

S  X  -  |Y  "  +1)  -  (2  +X)  ‘  (28) 

The  left-hand  equation  is  also  a  rewriting  of  equation  (13).  The  right-hand 
equation  does  not  involve  R.  This  means  that  the  two-parameter  family  of 
solution  curves  in  the  (R,  X,  Y)  space  have  a  one-parameter  family  of  projec¬ 
tions  in  the  (X,Y)  plane.  Among  the  latter  family  is  a  special  curve  which 
represents  the  constraint  between  Q  and  T  as  well  as  r  necessary  for  the 
pressure  to  become  vanishing  at  infinity.  Other  solutions  of  the  right-hand 
equation  are  unphysical,  for  the  subsequent  solutions  of  the  left-hand  equation 
have  the  defects  that  either  the  temperature  or  the  heat  flux  becomes  negative 
eventually.  On  the  other  hand,  the  left-hand  equation  is  linear  in  R.  This 
means  that  various  solutions  differ  only  in  scaling.  Ve  have  accounted  for 
the  scaling  by  the  inclusion  of  the  arbitrary  constant  (t  in  the  definition  of 
R  in  equation  (26). 
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It  turns  out  that  the  physically  meaningful  solutions  of  equations  (28) 
are  represented  by  three  curves,  which  emanate  from  the  critical  point  R^,  =  0, 
X^,  =  'T~  +  -^5  Yc  =  \  of  the  differential  equation  (28)  and 

recede  to  R  =  ».  The  one  with  X  >  X^,  ¥  >  Y^  has  the  asymptotic  property  that 
T  diminishes  like  r"2/7  for  large  r  (Parker,  1964).  The  one  with  X  =  X^,, 

Y  =  Yj,  has  T  diminishing  like  r"2/s  (Whang  and  Chang,  1965).  The  one  with 
X  <  Xj,,  Y  <  Yj,  has  T  diminishing  like  r"4/a  (Durney,  1971).  These  three 
curves  will  be  collectively  denoted  by  X  =  h(R) ,  Y  =  k(R) .  Their  other  projec¬ 
tions  will  be  denoted  by  Y  =  f(X)  and  R  =  g(X).  Clearly,  h(R)  and  g(X)  are 
inverse  functions  to  each  other,  whereas  k(R)  is  the  composite  function  of  f(X) 
by  h(R) .  Each  of  these  four  functions  has  three  branches  (see  Figure  1). 


The  two  branches  approach  X  =  X  ,  Y  =  Y  with  the  slope 


df 

dX 


x=x 


1  +  5A+  (2801  +  3410  X  +  1025  X2)1^ 
175  +  125  a. 


The  function  g(X)  has  the  following  power  series  representations. 


L- 3/ (7+5  A) 


ij(X)  =  <  indefinite 


1+  ro  ro<x,xc 


£e" 


for  X  =  X 


(30) 


1  + 


rr 

£  » 


for  X  <  ^ < 


The  coefficients  are  given  by  the  recurrence  formulas  (with  gQ  =  1,  gQ  =  1) 


15  n-1  >  7  i  _ 

’n  (14  +  10>)n  -^0*1  1  7  +  '  j;/  n-i  : 

n-1 

+  =  TZ  £  .  IV+  +  . 

ti  -tn  1=0  ;  2  /  n-i  i 


I  i  OC 

The  two  branches  approach  R  =  0,  X  =  X^,  in  a  manner  like  |X  -  Xj  with  the 
exponent 


,  5  df 

1  -  -  -rr-l 

dAt  X=X„ 


being  greater  than  unity. 


The  function  h(R)  has  the  following  power  series  representations 
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(  r-(7+5a)/3 

h(R)  =  X„ 


1+  »^-“l7+5S)/3 

n=l 


V 


,2/7 


i  *(i  ♦  f k/7-  E'2n/7 


for  X<  X 


for  X  =  X 


for  X  >  X. 


(31) 


with  the  coefficients  given  by  the  recurrence  formula  (with  ho  =  1,  hQ  =  1, 

+  _  35  245  jv 

hl  “  3S  +  “35"  K)  ‘ 

n-l 

2  y  V-  .  10+10A\(2n+3i)(7+5x)  +4l  +25XU. 

h»  -  -  sfesv1  +ii7r) - 45 - 1  “-1-1 


1 


n(2n+7) 


n-l 


.*  1 


J(6n-6l-35X)  \_±  +  £^(2n+3i)(2i+5)  h*  h+_.j  . 

i=i  J 


The  function  k(R)  has  the  following  power  series  representations: 


p  oo 


■n(7+5k)/3 


4(R)  =  <j  Yc 


for  y<  y 


for  Y  =  Y  (32) 

c 


i"2/7i  ?  -  Sk"R’2n/7i 


for  Y  >  y 


•  H 

with  the  coefficients  given  by  the  recurrence  formulas  (with  k  =  3  + 

kl=T 


v~  =  !  n(7  +  5  x  )  +  (10  +  5  X 


)  h' 

7 .  r 


<  "  i  <  ■”  +  5  >  hn 


•>Jro 


Now,  in  view  of  equation  (27),  the  compatible  heat  flux  is  given  by 


C  T5/2 

Q  -  K  p  u  T  f(  cos-x)  , 


(33) 


and  in  view  of  equation  (25),  the  additional  invariant  is  given  by 


1  c  t5/2 

7  K  rfu“  cos  '  }  • 


(34) 


With  these  results,  the  admissible  solution  can  be  calculated  in  terms  of  r 
PQ,  uq,  T  Bo,  and  XQ  in  the  following  algorithm.  First,  determine  the 
compatible  heat  flux: 


X 

o 

Y 

o 


i*  o  u 
O  i  o  o 


f(xj  , 

O 


cos 


% 

o 


y 


0. 

o 


=  K  Pc 


u  T  Y 

0  0  0  0 


Next,  evaluate  the  invariants: 

o 

7TI  =  - 


Do  uo 


B, 


£  =  ?ono(  1  2  _  ± 
^  Bn  \  2U°  2 


2  1  n  2  2  •  2  GM 

(1  ...  on  «  A  I 


ATr  sin^e  - 


k  +  H)  +  r  • 


=  ^ro  sin  ©j 


Po> 

B 

L  O 


(Uo  sin  \  +  %  r0  Sin  9)  -  A>  Sin^ 


7  =  r  “  B 


COS 


C  =  7  *(x  )  • 


Then  calculate  the  physical  quantities  at  the  heliocentric  distance  r  of 
interest : 
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R  =  Cr 


k  jn  7 _ h(R)  1 2//5 

C  cos^XQ  r  J 


1/2 

u  =  2  ' 


£  1,2  2 

—  +  —  r 

m  2  ° 


.  2. 
sxn  8  + 


GM 


-  |KT  -  KT  k(R) 
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/  1’  \ 
arctan,  —  tan  \  ) 
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B  = 
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r-  cos  A. 
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5.  Two-point  Relationship  in  Corotating  Streams 


Now  we  consider  a  corotating  solar  wind' stream  in  the  equatorial  plane 
(viz,  0  =  90°).  In  the  rest  frame,  this  spiraling  flux  tube  rotates  at  the 
angular  velocity  Q  =  360°/27  days.  Let  t  be  the  instant  when  this  corotating 

0  V 

stream  sweeps  through  Venus.  At  that  time,  the  earth  is  aziinuthally  ahead  of 
Venus  by  an  azimuthal  angle  given  by  4^.  -  We  want  to  know  the  later 

instant  t^.  when  the  same  stream  will  sweep  through  the  earth.  The  cross- 
section  (denoted  by  E'  in  Figure  2)  of  the  corotating  stream  that  sweeps 
through  the  earth  is  slightly  azimuthally  behind  the  cross-section  that  sweeps 
through  Venus.  The  azimuthal  angle  between  these  two  cross-sections,  one  has 
the  heliocentric  distance  r^  =  0.72  AU  and  the  other  r^.  =  1.0  AU  is 

tan(-V) 

%  '  V  =  ^rE  '  V  * 


in  which  Xy  is  the  spiral  angle  of  the  flux  tube  at  the  heliocentric  distance 
ry  Equation  (35)  is  obtained  from  the  differential  equation  rdij>/dr  =  B^/B^ 
for  the  field  line  together  with  the  equation  B^/B^.  cs  (r/r^)  tan  Xy  along  the 
flux  tube.  Since  the  earth  is  revolving  aroun  the  sun  at  the  angular  velocity 
ftj,  =  360°/ 365  days,  the  time  delay  tg,  -  ty  is  equal  to  (<j>g.  -<J -fig.). 

Hence 


t 


E 


_E 

i. 

o 


tan(-iy) 


(36) 


Note  that  Xy  is  negative,  for  angles  are  reckoned  positive  in  the  eastward 
direction  of  the  solar  rotation. 
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Next,  suppose  that  the  observed  solar  wind  data  at  Venus  (as  measured  by 

the  space  probe  Pioneer-Venus)  have  the  mass  density  Py,  the  bulk  speed  Uy, 

the  temperature  Ty,  the  magnetic  strength  By  and  the  spiral  angle  Xy  From 

the  transformation  =  u  cos  x»  =  u  sin  x  +  between  the  flow  velocity 
— ► 

U  in  the  rest  frame  and  the  flow  velocity  u  in  the  corotating  frame,  the  flow 
speed  in  the  corotating  frame  is  given  by 


w 


A)1'' 


+  ry  sin(-\) 


(37) 


Using  the  algorithm  described  in  the  previous  section,  we  can  calculate  the 
mass  density  p£, ,  the  flow  speed  Ug,,  the  temperature  Tg, ,  the  magnetic 
strength  Bg,  and  the  spiral  angle  Xg.  at  the  heliocentric  distance  Tg.  The 
corresponding  bulk  speed  in  the  rest  frame  is 


UE-  =  (V  +  "E  V  aiU  \+  V  rE  f' "  ■  (38) 

These  extrapolation  values  are  to  be  compared  with  the  observation  values  of 
the  solar  wind  at  the  earth  at  the  instant  t£  given  by  equation  (36).  The 
results  will  indicate  to  what  extent  the  apparent  temporal  variation  in  the 
interplanetary  space  is  due  to  the  effect  of  corotation.  From  such  a  study, 
we  can  infer  the  time  scale  in  which  the  temporal  effect  due  to  the  intrinsic 
short-period  changes  of  coronal  condition  may  be  neglected.  In  the  following  se 
tion  we  shall  present  the  results  obtained  from  the  comparison  of  the  Pioneer- 
Venus  data  with  the  data  of  earth-bound  satellite  ,  based  on  the  model  described 


in  this  section. 


Correlation  Features  in  Observational  Data 

We  use  observational  data  of  Pioneer  12  and  IMP  8  for  the  period  from 
January  1980  to  June  1980  to  examine  the  two-point  correlation  features  due 
to  corotation  of  the  solar  wind  streams  in  the  interplanetary  space.  During 
this  period,  Venus  to  which  Pioneer  12  is  bound  was  azimuthally  behind  the 
Earth  to  which  IMP  8  is  bound.  The  Earth-Sun-Venus  angle  decreased  from  103° 
on  January  1,  1980  to  0°  on  June  15,  1980.  Subsequently,  Venus  became 
azimuthally  ahead  of  the  Earth,  by  120°  on  December  31,  1980.  An  azimuthal 
separation  of  lr  amounts  to  one  day  for  the  delay  time  of  the  rotating  stream 
tubes . 

Typical  values  of  solar  wind  speed,  density  and  temperature  observed 
by  Pioneer  12  were  published  in  digital  format  in  Solar-Geophysical  Data, 

Part  I  (N0AA,  1980),  whereas  variations  of  solar  wind  speed,  density  and 
temperature  observed  by  IMP  8  were  published  in  graphic  format  in  Solar- 
Geophysical  Data,  Part  II  (N0AA,  1980).  Each  triplet  of  data  from  Pioneer  12 
represents  the  property  of  a  solar  wind  stream  at  0.7  AU.  We  use  equations 
in  section  ♦,  ,ith  /.*=  .  ,to  calculate  the  values  of  speed,  density  and 
temperature  at  1.0  AU.  We  also  calculate  the  delay  time,  using  Equation  (Hi). 
The  extrapolated  values  of  speed,  density  and  temperature  at  delayed  times 
were  plotted  as  dots  in  Fi.-u:-  ?.  The  segmental  curves  in  the 

figures  display  the  observational  data  of  IMP  8.  Coincidence  of  the  three 
dots  representing  the  extrapolated  triplet  with  the  three  curves  would  mean 
perfect  correlation.  Since  the  data  by  IMP  8  were  not  in  digital  format, 
the  correlation  feature  is  to  be  examined  visually  rather  than  to  be  assessed 
in  quantitative  manner. 


7.  Discussions  and  Conclusions 


In  this  study  we  have  examined  the  data  from  Pioneer  12  and  IMP  8  to 
discern  the  corotational  correlations  between  the  two  sets  of  data  observed 
at  different  regions  of  the  interplanetary  space.  A  high  correlation  is 
expected  when  the  solar  wind  is  in  a  quasi-steady  state.  We  compare 
theoretically  extrapolated  values  from  Pioneer  12  data  with  corresponding 
IMP  8  data.  The  correspondence  is  determined  by  the  assumption  that  solar 
wind  streams  corotate  with  the  Sun.  The  extrapolation  is  done  through  a 
suitable  modeling  for  solar  wind  streams  in  the  far  region  of  the  inter¬ 
planetary  space. 

The  data  from  Pioneer  12  and  IMP  8  as  published  in  Solar-Geophysical 
Data  (NOAA,  1980  )  consist  of  the  three  values  of  the  flow  velocity, 

particle  density  and  temperature  for  the  proton  species,  together  with  the 
azimuthal  separation  angle  between  the  Earth  and  Venus,  at  selected  times. 

We  use  the  standard  one-fluid  solar  wind  equations  with  thermal  conduction 
to  extrapolate  the  velocity,  density  and  temperature  at  1.0  AU  from  the 
observed  values  at  0.7  AU. 

The  model  we  have  used  for  solar  wind  streams  is  capable  of  accounting 
for  correlation  features  arising  from  corotation  of  solar  wind  streams  alone. 
There  are  many  physical  factors  in  the  model  whose  simple-minded  treatment 
cause  deviations  from  good  correlation  between  theoretically  derived 
quantities  and  observationally  obtained  quantities.  Among  them  are  non¬ 
steady  coronal  conditions  and  non-radial  streamlines.  Each  of  these  two 
is  beyond  the  present  state  of  arts  for  remedy.  This  is  so  not  only  because 
of  the  mathematical  complexities  in  any  time-varying  or  global  treatment, 
but  also  because  of  the  observational  lack  of  structured  coronal  conditions. 

Generally  speaking,  better  correlation  between  interplanetary 


points  is  achieved  when  the  azimuthal  separation  between  them  is  small. 

This  is  expected,  for  the  simple  reason  that  less  temporal  changes  in  the 
coronal  condition  enter  the  solar  wind  streams  in  a  shorter  duration  of 
delay  time. 

In  the  realistic  interplanetary  space,  temporal  changes  are  always 
present.  It  seems  plausible  to  regard  these  temporal  changes  as  a  series 
of  quasi-steady  states.  If  so,  the  relevant  question  is:  what  is  the  time 
scale  over  which  the  averages  behave  like  quasi-steady  quantities.  For 
the  purpose  of  a  prediction  scheme,  the  time  scale  should  increase  as  the 
azimuthal  separation  angle  increases.  When  the  time  scale  exceeds  the 
characteristic  time  of  the  magnetospheric  phenomena  of  interest,  predic¬ 
tion  no  longer  makes  sense.  To  assess  the  time  scale  of  acceptable 
correlation,  it  seems  worthwhile  to  perform  correlation  study  based  on 
temporally  averaged,  with  chosen  time  scales,  data.  The  present  study  uses 
instantaneous  values  of  the  interplanetary  conditions  at  Venus  and  at  the 
Earth.  With  the  efficient  methodology  and  numerical  code  developed  in  this 
study,  the  suggested  study  should  be  an  accomplishable  task. 
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Figure  Captions 

Figure  1.  Plots  of  the  functions  f(X)  and  g(X).  They  describe  the  (Y,X) 
projection  and  the  (R,X)  projection  of  the  critical  solution  of  equa¬ 
tion  (28)  with  A  =  0. 

Figure  2.  A  corotating  solar  wind  stream  in  the  equatorial  plane.  The 
stream  flows  along  a  flux  tube  emanating  from  the  sun  (denoted  by  S). 

It  sweeps  through  Venus  (denoted  by  V)  at  time  ty  and  through  the  earth 
(denoted  by  E)  at  a  later  time  t^.  The  flux  tube  has  a  spiral  angle  Xy 
at  the  heliocentric  distance  r^  of  Venus.  The  cross-section  (denoted  by 
E')  sweeping  through  the  earth  has  the  heliocentric  distance  r^  of  the 
eartu  and  is  behind  the  cross-section  sweeping  through  Venus  by  an  azi¬ 
muthal  angle  <}»y  -  $  ,  .  The  earth  is  ahead  of  Venus  by  an  azimuthal  angle 

4>E  "  3t  time  t y. 

Figure3.  Time-sequence  plots  of  the  interplanetary  conditions,  indicated 
by  particle  density,  flow  speed  and  temperature,  of  the  Earth.  Ti.?  dots 
depict  the  theoretical  values  extraplolated  from  the  observational  data  of 
Pioneer  12,  whereas  the  curves  depict  the  observational  data  of  IMP  8. 

The  Earth-Sun-Venus  angle  is  indicated  on  the  top  of  panel. 
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